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Abstract 

In this paper statistical solutions of the 3D Navier-Stokes-a model with periodic boundary condition are 
considered. It is proved that under certain natural conditions statistical solutions of the 3D Navier-Stokes-a 
model converge to statistical solutions of the exact 3D Navier-Stokes equations as a goes to zero. The 
statistical solutions considered here arise as families of time-projections of measures on suitable trajectory 
spaces. 
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1 Introduction 

In this work we investigate the convergence of statistical solutions of the three-dimensional Navier-Stokes-a 
model to statistical solutions of the three-dimensional Navier-Stokes equations, as a goes to zero. We con- 
sider the equations with periodic boundary conditions and zero space average, and the statistical solutions are 
considered in the sense recently introduced by Foias, Rosa and Temam in [25l [23] . 

Most of the knowledge concerning turbulent flows are rooted in heuristic and phenomenological arguments. 
One of the fundamental observations is that, although quite irregular, turbulent flows display a certain order in a 
statistical sense, so that mean quantities and other low order moments are usually more regular and predictable. 
The statistical solutions that we investigate are mathematical objects used to address the statistical properties 
of the flow in a rigorous mathematical way directly from the equations of motion. 

The Navier-Stokes equations have been widely used as a model for Newtonian turbulent flows (see e.g. 
[2] [551 ISH 13H] ) and a rigorous mathematical formalization of the statistical theory of turbulence came with the 
introduction of statistical solutions for these equations. Several statistical estimates can be obtained rigorously 
using statistical solutions; see for instance [16] [TO] [3] [17] [27J [20j [34]. Therefore, a better understanding of 
statistical solutions is crucial for a rigorous mathematical theory of turbulence. 

The concept of statistical solutions was first introduced by Foias and Prodi [13 , 22 in the early 1970's (see also 
an earlier related mathematical work by Hopf [3UJ ) . They considered as a statistical solution a family of measures 
on the phase space satisfying a Liouville-type equation together with some regularity conditions. Some years 
later Vishik and Fursikov [43] introduced a different type of statistical solutions, given by measures on suitable 
trajectory spaces (see also [UJES]). More recently, in [2SJ[23], Foias, Rosa and Temam elaborated a notion of 
statistical solutions which, in some sense, relates the two previous notions and has better analytical properties 
than the previous ones. More precisely, they considered a Borel probability measure on the space C([0, T], H w ) 
which is carried by the trajectory space of Leray-Hopf weak solutions. The family of time projections of such 
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a measure gives rise to a statistical solution in the sense of Foias-Prodi that possess additional analytical 
properties. 

We also mention a work due to Capinski and Cutland, [5], preceding the works of Foias, Rosa and Temam, 
in which the authors prove, using non-standard analysis, the existence of space-time statistical solutions (and 
of individual weak solutions) defined in the same spirit as that in |25[ 123] , with the main difference that they 
use a slightly different definition of weak solutions upon which the statistical solutions are built. 

There is also a large literature on the stochastic version of the Navier-Stokes equations, which relics on some 
ideas from the Vishik-Fursikov formulation; see e.g. |13j . 

A few other equations have also been considered as suitable models for turbulent fluid motions in specific 
aspects. For instance, the Navier-Stokes- a model is known as a good approximation model for well-developed 
turbulent flows (see [HI |5] for the cases of infinite pipes and channels). Also, the 3D Navier-Stokes- Voigt 
equations are known as an appropriated model for direct numerical simulations of turbulent flows in statistical 
equilibrium. For a survey on approximating models of the Navier-Stokes equations and their properties we refer 
the reader to the Introduction given in [55] . Regarding the importance of statistical solutions to yield statistical 
estimates rigorously, we highlight the work of Ramos and Titi, |35j . The authors derived statistical properties 
of the invariant measures associated with the solutions of the 3D Navier-Stokes- Voigt equations and established 
the convergence of probability invariant measures associated with the 3D Navier-Stokes- Voigt equations to 
stationary statistical solution of the 3D Navier-Stokes equations. With this result, they argue, via statistical 
estimates obtained rigorously, that the 3D Navier-Stokes- Voigt is in fact a reliable subgrid scale model for direct 
numerical simulations of turbulent flows. 

We focus our study on the Navier-Stokes-a model and on the statistical solution defined by Foias, Rosa and 
Temam, namely a family of time projections of a measure carried by the trajectory space of the Leray-Hopf 
weak solutions; see Section 12.61 The Navicr-Stokcs-a model (also known as Camassa-Holm equations) was 
introduced by Chen et al, in [B]. This model is a regularized approximation of the 3D Navier-Stokes equations 
such that, in some terms of the equation, the velocity field is replaced by a smoother (filtered) velocity field 
depending on a small parameter a > 0; see Section 12.41 As observed in |I8j . this regularized approximation 
introduces an energy penalty that inhibits the creation of smaller and smaller excitations below the length scale 
a. More precisely, it was proved in [T5] that the wavenumber spectrum of the translational kinetic energy for 
the Navier-Stokes-a model rolls off as k~ 3 for na > I instead of continuing along the Kolmogorov scaling law, 
k~ 5 / 3 , which is followed for na < I. 

In this article we prove that, under certain natural conditions, statistical solutions of the Navier-Stokes-a 
model converge to statistical solutions of the Navier-Stokes equations as a goes to zero. We consider both 
stationary and time-dependent statistical solutions. We point out the importance of stationary statistical 
solutions in the study of turbulence in statistical equilibrium in time, while, time-dependent statistical solutions 
are useful in the study of evolving or decaying turbulence. 

Another important fact is that the notion of stationary statistical solutions provides a generalization of the 
notion of invariant measures of semigroups. In the case of the 3D-Navier-Stokes equations this is currently 
needed since the well-posedness of the equations has not been established. In the two-dimensional case, in 
which the Navier-Stokes equations have a well-defined semigroup, the two notions have been proved to agree 
with each other. 

The definition of statistical solution for the Navier-Stokes-a model is inspired by the corresponding definition 
for the Navier-Stokes equations, being the family of projections in time of a Borcl probability measure in a 
suitable trajectory space and carried by the set of individual solutions of the equation. The natural trajectory 
space for the Navicr-Stokcs-a model is C([0, T], H), which is included in C([0, T], H w ), which is the natural space 
for the Navier-Stokes equations, so that both statistical solutions can be regarded as projections of measures on 
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the same space. In this way, both statistical solutions can be viewed as mathematical objects of the same type, 
allowing a direct comparison between them, and yielding a natural framework for studying the convergence of 
the corresponding statistical solutions as a goes to zero. Of course, since the Navier-Stokes-a model is well- 
posed, its solution semigroup induces a measure in the trajectory space, starting from any initial measure in 
the phase space, in such a way that the time-projections of such a measure form indeed a statistical solution 
in the sense we consider here (see Section 2.6). It is expected that the converse is also true, namely, that 
any statistical solution for the Navier-Stokes-a model is the family of projections of a measure induced by the 
solution semigroup (as it happens for the two-dimensional Navicr-Stokcs equations), but we do not address this 
issue here. 

This is one of a few results about convergence of statistical solutions (see e.g. [T2J[TTJ[35] for other equations) 
and it is the first result of convergence for this type of statistical solution, and we believe the main ideas presented 
here will be valuable for extending the result to other types of approximations (see [4]). 

The structure of the paper is as follows. We start with some usual definitions and basic results concerning 
the functional-analytic framework. In Section 12.21 we review some facts regarding Borel measures and we state 
some compactness results in the space of Borel measures which may not be so familiar to the reader and which 
will be used to obtain the convergence of the measures related with the statistical solutions. Later, in Sections 
12.31 and 12.41 we briefly introduce the Navier-Stokes equations and the Navier-Stokcs-a model, focusing on the 
results that will be important throughout this paper. In Sections l2.5l and l2"i)I we recall the definition of statistical 
solutions for the Navier-Stokes equations and introduce a notion of statistical solution for the Navier-Stokes-a 
model. 

The main goals of this paper are presented in the final Section [3l where we state and prove the convergence 
theorems for the measures on the trajectory space (Theorem 13.11 and Corollary I3.3[) and for the statistical 
solutions fThcorem l3.2[ Corollary |3T4"l and Theorem I3.3|) . This section is divided into two parts, one for the time- 
dependent case and the other for the particular case of stationary statistical solutions. The main condition for 
the convergence of a family of statistical solutions of the a- Navier-Stokes equations to a statistical solution of the 
Navicr-Stokcs equations as a goes to zero is that the mean kinetic energy of the family be uniformly bounded. 
This provides the tightness of the family of measures that guarantees the compactness of the family according 
to the theory of Topsoe discussed in Section 12.21 The main work, then, is to show that this compactness is 
sufficient to guarantee that the limit family of measures is a statistical solution of the Navier-Stokes equations. 
A crucial step in this proof is to show that the limit measure is carried by the space of Leray-Hopf weak solutions 
of the Navier-Stokes equations. This step depends very much on a result by Vishik, Titi, and Chepyzhov |46] 
on the convergence of individual solutions. 

2 Preliminaries 

In this section, we set the notation and provide the definitions and results that are needed throughout this work. 
2.1 Functional setting 

Let Q, := llf =1 (0,ii), where L t > 0, for i = 1,2,3, and let C~ r (rj;M 3 ) represent the space of the infinitely 
diffcrentiable functions u, with u(x) = (iti(x), it2(x), «3(x)) £ M 3 defined for x = (x\,X2,X3) £ M 3 , which are 
f2-periodic. We define the set of periodic test functions with vanishing average as 
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the norm in H are defined, respectively, by 

(u,v) := / u • vdx and |u| := \J (u, u), 
Jn 

where u • v = u»Ut, and in V they are defined by 

((u,v)) := (Vu,Vv) and ||u|| := y/((u,u)), 

where it is understood that Vu = (dui/dxj) 3 - =1 and that the second term is the integral of the componentwise 
product between Vu and Vv. We also consider the space H endowed with its weak topology, which we denote 
by H w . 

Let A be the Stokes operator defined as A = —PA, where P : L 2 (fl) 3 — > H is the Leray-Hclmholtz 
projection, i.e., the orthogonal projector in L 2 (f2) 3 onto the subspace of divergence-free vector fields. We 
denote by D(A) the domain of A, which is defined as the set of functions u G V such that Au G H. Recall 
that, in the periodic case with zero space average, ^4u = — Au for u G D(A) = V n 2 (SI) 3 and A is a 
positive self-adjoint operator with compact inverse, so that it has a sequence {Ai}i 6 N of positive eigenvalues 
counted according to their multiplicity, in increasing order, associated with an orthonormal basis {w^jigN in 
H. Furthermore, the Poincare inequality holds, i.e., for all u G V, 

Ai|u| 2 <||u|| 2 , (1) 

where Ai > is the first eigenvalue of the Stokes operator. 

Denote by P k the Galerkin projector defined as the projector onto the space spanned by the eigenfunctions 
associated with the first k eigenvalues, i.e., 

k 

P k u = y^(u, w,)w„ Vu 6 H. 

i=l 

There are some important properties of P k that we want to highlight. For instance, we have that, for all u G H 
and for all k G N, |Pfcii| < |Pfc+iu| and |Pfcu| < |u|. We also have that P& : H w — > H is continuous. Indeed, 
observe that the open sets in H can be characterized by a basis of neighborhoods 

C(u, r) = {w e H : |u - w| < r}, 

for u G H and r > 0. On the other hand, the collection of open sets in H w has a characterization by a basis of 
neighborhoods given by 

N 

0„(u,r,vi,...,vjv) = {w G H w : ^|(u-w,Vj)| 2 < r 2 }, 

i=l 

for u G H w . r > 0, TV G N, and vi, . . . , vjv G H. Notice that if we prove that for all u G H w and r > the 
set P k ~ 1 (0(PkU, r)) is open in H w , then we obtain that Pk : H w — > H is continuous. This follows from the fact 
that P,7 1 (0(PkU, r)) is of the form 

P^{0{P k u,r)) = {w G H w : P fe w G 0(P k u,r)} = {w G H w : \P k u - P k w\ < r} 

k 

= {wGP tu :^|(u-w,w J )| 2 <r 2 }, 

»=i 

and, hence, is an element of the basis for H w . 
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For the functional setting concerning the Navier-Stokes-a model, we adopt the framework introduced by 
Vishik, Titi and Chepyzhov in 

We start with the natural space for the solutions of the Navier-Stokes-a model. Given an interval Jcl, 
we define 

Tx = {z : z(-) e Ll c (I; V) n L% C (I; H), d t z(-) E L 2 oc (I; D(A)% (2) 

We endow this space with its natural weak-type topology, which we term r topology, and which can be defined 
in terms of nets as follows: a net of functions {z 7 } 7 C J 7 / converges to a function z G Ti in the topology r if 
for each compact interval J d I, 

z 7 A z in L°°{J;H), z 7 z in L 2 ( J; V), and <9 f z 7 — ^ d t z in L 2 (J; D(A) 1 ). 

Consider also the following Banach space 

T\ = {z : z(-) e L\(I- V) n L°°(J; H),d t z(-) G L 2 h (I; D(A)')}, (3) 

with norm given by 

ll z ll^b = II z IIl^(/,v) + \\z\\l<~(i,h) + \\dtz\\ L 2 {ID{Ay) , 

where 

i-t+i 

Ml?(i.v)= sup / ||z(s)|| 2 ds, 
{te/:t+ie/} Jt 

and 

/•t+i 

\\9M\lI(i,d(AY) = SU P / ll 9 i z ( s )lll>(A)' rfs - 
{tei-.t+iei} Jt 

Finally, we introduce a natural space for the solutions of the Navier-Stokes equations, which is the space 
Ci oc (I,H w ) of continuous functions from an interval I C R to H Wl where, as defined above, H w stands for 
the space H endowed with the weak topology. This space can also be seen as the space of weakly continuous 
function from / to H. The topology on Ci oc (I , H w ) is that of uniform convergence in H w on compact intervals 
oil. 

2.2 Measures and Compactness on the space of measures 

The natural measure space in our framework is the space of Borel probability measures on Ci oc (I , H w ^j. Since 
we are interested in the study of the convergence of family of measures we need a compactness result on the 
space of measures. For measures on a separable metrizable space there is the well-known compactness result 
due to Prohorov [33]. In our case, however, Ci oc {I, H w ) is not metrizable, so we use a more general compactness 
result due to Topsoe (see [101 HB H2), which is suitable for our framework. 

Let X be a Hausdorff space and let Bx be the Borel a- algebra on X. The set of finite Borel measures on X 
is denoted by M(X) and the set of Borel probability measures is denoted by V{X). We say that a measure /i 
is tight if for every set A £ Bx, 

n(A) = sup{n(K) : K is compact, K C A}. 

Denote by M.{X; t) the set of all Borel finite measures which are tight, and by V(X; t) the set of all measures 
in V(X) which are tight. If X is a Polish space (i.e., a separable and completely metrizable space) then every 
finite Borel measure is tight, so that M(X;t) = M{X) and V(X;t) = V(X). 

We say that a net (/x 7 ) in V(X) is uniformly tight if for every e > there exists a compact set K C X 
such that 

t x 1 {X\K)<e, V 7 . 
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Let Y be a Hausdorff space and let F : X —> Y be a continuous function. For a measure /j, G A4(X) wc 
define the measure F/j, induced by fi on Y as Ffi(E) = ^{F^ 1 {E)), for every Borel set E C Y. In this setting, 
the Change of Variables Theorem (see e.g. [1]) says that if tp : Y — > M is a F/U-integrable function then 93 o F is 
^t-integrable and 

/ <p(F(x))dn(x) = f V (y)dFfi(y). (4) 
Jx Jy 

Now, for the result of compactness on the space of measures we endow the space M(X) with the weakest 
topology for which the mapping \i i-> fJ,(f) is upper semicontinuous for every / bounded, real-valued, upper 
semicontinuous function on X, where n(f) stands for the integral f x f(x)d/i(x). We use the symbol — » to 
denote the convergence of nets in M.(X) with respect to this weak topology. The spaces V(X), A4(X;t), and 
V(X; t) are endowed with the topology inherited from Ai(X). 

Recall that a topological space X is completely regular if every nonempty closed set and every singleton 
disjoint from it can be separated by a continuous function. 

In [50], Topsoe proved a result of compactness on the space of measures on an abstract space X. In the case 
X is a topological Haudorff space the result is reduced to the following(see |41[ Theorem 9.1]): 

Theorem 2.1. Let X be a Hausdorff space. Let (/i 7 ) be a net in V(X\f) which is uniformly tight. Then, there 
exist ji G V(X;t) and a subnet (// 7/3 ) such that fi^ p A fi. 

Topsoe also proved, in [5T], the following result: 

Lemma 2.1. Let X be a completely regular Hausdorff space. For a net (/x 7 ) in M.(X) and /i G M(X,t), the 
following statements are equivalent 

(1) n-i -> 

(2) limsup^ 7 (/) < fJ,(f), for all f bounded upper semicontinuous function; 

(3) liminf /x 7 (/) > fJ,(f), for all f bounded lower semicontinuous function; 

(4) lim/z 7 (/) = fJ-(f), for all bounded continuous function f. 

Remark 2.1. The last lemma was actually stated in a more general setting. For instance, if X is only Hausdorff 
and \i G M(X) then the first three statements are equivalent and each of them implies the last one. If X is a 
completely regular Hausdorff space and fj, G M.(X) is r -smooth (a condition which is satisfied for every tight 
measure) then all the statements are equivalent. 

Observe that, by Theorem 12. 1[ if (/i 7 ) is a net in V(X; t) which is uniformly tight then there exists a subnet 
that converges to a limit \i which is tight. Moreover, by Lemma [2.11 if X is a completely regular Hausdorff 
space, the convergence in the weak topology in V{X]t) is equivalent to the usual convergence // 7(3 (/) — > ^{f), 
for every / G Cb(X), denoted by /i 7/3 [i, where Cb(X) denotes the set of bounded, continuous, real-valued 
functions defined on X. 

We state the last observation as the next theorem: 

Theorem 2.2. Let X be a completely regular Hausdorff space. Let (/i 7 ) be a net in V(X;t) which is uniformly 
tight. Then, there exist [i G V(X;t) and a subnet (/i 7/3 ) such that /x 7/3 — \i, i.e., 

lim M7i3 (/) = M/)> for all f eC b (X). 
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Also in [5T], Topsoe stated the following result, which we prove here with more details. 
Theorem 2.3. Let X be a H aus dor ff space. Then, V(X;t) is a Hausdorff space. 

Proof. First, recall that a Hausdorff space can be characterized as a topological space where every net converges 
to at most one point. Therefore, it is enough to prove that if (^ 7 ) 7 is a net in V(X;t) which converges to two 
measures fii,/i2 £ V(X;t), i.e., fi 7 — > /i\ and /i 7 A /i2, then /i\ = /i2- Let A £ Bx, denote by A the interior 
of A and by A the closure of A. It is clear that the characteristic functions \a anc ^ %A are, respectively, lower 
semicontinuous and upper semicontinuous functions. Therefore, using Lemma l2.1[ we obtain that 

Hi(A) < liminf fa (A) < limsup/z 7 (^4) < /12(A). 

7 7 

Now, let E S Bx, and let us prove that /i\(E) < /12(E). In order to do so, consider any compact sets K\ C E 
and K2 C E c . Since X is Hausdorff there exist disjoint open sets A and B such that such that K\ C A and 
A 2 C 5. It is clear that A C X \K 2 . Thus, 

Mi(A'i) < fix(A) < /12(A) < /i 2 (X \K 2 ) = l- /i 2 (K 2 ), 

which leads us to 

/i 1 (A 1 )+/i 2 (A 2 ) < 1. 

Since K\ and K 2 are arbitrary compact sets satisfying K± C E and K2 C E c , we can take the supremum over 
all compact sets K± C E and the supremum over all compact sets K2 C E c in the last expression, and we find 
that 

sup{/ii(Ai) : K\ is compact, K\ C A} + sup{^ti(A2) : K2 is compact,^ C E c } < 1. 
Since /i\ and ^2 are tight, we conclude that 

/H(E) + /i 2 (E c ) < 1. 

Thus fii(E) < /12(E), for all E £ Bx- Now, since /i\(X) = /i2(X) = 1 it follows that /i\ = /i 2 . □ 
Let A" be a completely regular Hausdorff space and /i\,/i2 & V(X;t). Then, 

/ii = /12 if and only if / (p(x)d/ii(x) = / ip(x)d/i2(x) 7 Vip £ Cb(X). (5) 
J x Jx 

Indeed, suppose that J x tp(x)d/i\(x) = J x (p(x)d/i2(x), for all ip £ Cb(X). We define the net (jU 7 ) 7 , where 
/i 7 = /i2, for all 7. Then, it is clear that /i 7 /i\ and /i 7 /i2- Since P(X;t) is Hausdorff we find that 
A*i = fi>2- The other implication is trivial. 

2.3 The Navier-Stokes equations 

We state only the results, properties and estimates that are needed throughout this work. For a more complete 
theory of the Navier-Stokes equations the reader is referred to [HI [HI I3U I3Z1 13H] , and the references therein. 

In the estimates below, we will consider certain quantities which are constant with respect to the solutions 
of the equations, but which may depend on the coefficients of the equations, the forcing terms and the spatial 
domain. We will call them non-dimensional constants when they are independent of re-scalings of the equations 
in space and time, hence, in particular, they may depend on the shape of the domain, but not on the size of 
the domain. Moreover, a non-dimensional constant will be called universal when it does not depend on any of 
the parameters of the equations, not even on the shape of the domain. 
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We recall that the Navier-Stokes equations can be written in the following functional form 

u t + vAa + B(u, u) = f , (6) 

where A is the Stokes operator and B(u, u) = P[(u • V)u], for u 6 V. 

The notion of solution that is considered here is the well-known Leray-Hopf weak solution, which is defined 
below. 

Definition 2.1. Let f £ Lf oc (I,V'). A function u is called a Leray-Hopf weak solution if: 

(i) u G L% C (I, H) n Ll c (L, V) n C loc (L, H w ); 
(ii) d t ueL^(I,V>); 

(Hi) u satisfies the weak formulation of the Navier-Stokes equations, i.e., 

u t + vAw + B(u, u) = f , 
in V' , in the sense of distributions on I; 
(iv) u satisfies the energy inequality in the sense that for almost all t' E I and for all t £ I with t > t' , 

i|u(t)| 2 + ^Vwi| 2 d.s<i|u(0| 2 +^(f( S ),u( s ))d S ; (7) 

(v) If I is closed and bounded on the left, with left end point to, then the solution is strongly continuous in H 
at to from the right, i.e., u(t) — > u(to) in H as t — > t^ . 

The set of allowed times t' in (JT]) can be characterized as the points of strong continuity of u, in H , from 
the right. In particular, condition (v) implies that t 1 = to is allowed in that case. 

Suppose that f S L°°(I, H) and let u be a Leray-Hopf weak solution. It is known that (see e.g. [HI Appendix 
ILB.l]) 

|u(t)| 2 < |u(0| 2 e- Al ^- 4 ') + J^||f HW Wl ~ (8) 

for all t' £ I allowed in © and for all i € / with t > t' . 

Furthermore, for all t' £ I allowed in (0 and for all i E J with t > t' , 



/ /■« \ 1/2 

1 /2 

( f \\d t u(s)\\l iAy ds) < — £1— |u(t')| 2 + 4tJ M i + ^ 5/2 A! /4 A/i(< - t'), (10) 

where c± is a universal constant and Mi is a non-dimensional constant which depends only on non-dimensional 
combinations of the parameters v, X\ and ||f ,h)- 
Define 

R* ~ ^-Ml-^h) (11) 

and observe that if |u(t')| < R, for some R > Rq and some t' £ I allowed in (J7J, then it follows from © that 
|u(t)| < R for all t > t'. 



(9) 



and 
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2.4 The Navier-Stokes-a model 

Most of the definitions, properties and results described below were taken from |46| . 
We consider the 3D Navier-Stokes-a model, in the periodic domain f2: 

v t — vAv + u x (V x v) + \7p = f , 

v = u-a 2 Au, (12) 
V • u = 0, 

where u is the unknown (filtered) velocity field, v is an auxiliary variable, p is the pressure, f is the external 
force, and a > is a constant. 

A functional formulation of (fT2]) can be written as 

— (u + a 2 Au) + uA(u + a 2 An) + B(u, u + a 2 Au) = f, (13) 
dt 

where A is the Stokes operator, and B(u, v) := P[u x (V x v)] is defined for u,v £ V, with values in V 1 , and 
which can be extended continuously to an operator from (u, v) 6 V x H with values in D(A)' (see [T§]). 
The notion of solution for the Navier- Stokes- a model that we consider is stated in the next definition: 

Definition 2.2. Let f G L 2 oc (L, H). A function u is a solution of U2)) on I if 

(i) u€L% c (I;V)nLl c (I;D(A)); 

(ii) d t ueL 2 oc (I-H); 
(m) u e Ci oc (I;V); 
(iv) u satisfies 

^■(u + a 2 An) + vA{u + a 2 Au) + B(u, u + a 2 Au) = f 
in D(A)' , in the sense of distributions on I; 
(v) u satisfies the energy equality in the sense that for all t',t S I with t > t' , 

i(|u(<)| 2 + a 2 \\u(t)\\ 2 )+v A||u( S )|| 2 + a 2 \Au(s)\ 2 )ds 



(14) 



= i(|u(0| 2 + a 2 ||u(i')|| 2 ) + J (f(s),u(s))ds. 
Observe that conditions (ii), (Hi) and (v) are consequences of (i) and (iv). 

The existence and uniqueness theorem of solution to the Navicr-Stokes-a model was proved in |19j : 

Theorem 2.4. Let f S H and Uq G V. Then, for each T > 0, there exists a unique solution of 112(1 on [0,T] 
in the sense of Definition ] 2. ^1 with initial data Uq. 

It is not hard to see that, in fact, one can assume f £ L°°(I, H) and still have the same conclusion as in the 
previous theorem. 

Another important property, which was observed and used in [46] is the following: If u is a solution of the 
Navier-Stokes-a model on /, then w, defined by w = (1 + a 2 A) 1 / 2 u, satisfies the functional equation 

w t + vAvr + (1 + a 2 A)~^ 2 B((l + a 2 A)^ 1 / 2 w, (1 + a 2 A) 1 ' 2 ™) = (1 + a 2 A)~^ 2 f, (15) 

and vice-versa. 

Therefore, based on the definition of solution to the equation (fT2|) . a notion of solution for the equation (ITSl) 
can be defined as follows: 
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Definition 2.3. Let f G Lf oc (L;H). A function w is a solution of M5\) on I if: 

(i) v,£L% c (I;H)nLl c (I;V); 

(ii) d t w€Ll c (I;D(A)>); 
(Hi) w E Ci oc {I-,H); 

(iv) w satisfies 

w t + i^Aw + (1 + a 2 ,4)- 1 / 2 i?((l + a 2 /!)- 1 /^, (1 + a 2 A)^ 2 w) = (1 + a 2 A)- 1 / 2 f , 
in D(A)' , in the sense of distributions on I ; 
(v) w satisfies the energy equality in the sense that for all t',t G I with t > t' , 

i|w(i)| 2 + v J* \\w(s)\\ 2 ds = l\MtW + j\(l + * 2 A)- 1/2 f(s)Ms))ds. (16) 

Again, here we have that conditions (ii), (Hi) and (v) are consequences of (i) and (iv). 

Suppose that f G L°°(I, H). We observe that if w is a solution of the Navier-Stokes-a model on an interval /, 
in the sense of Definition 12. 31 then |w(-)| 2 is an absolutely continuous function on / (see [461 Corollary 2.1]). As 
a consequence, one can show that, for any -0 : [0, oo) — > R, such that ip G C 1 ([0, oc)), xjj > 0, sup r>0 ip'(r) < oo, 
w satisfies the following estimate, 

V;(|w(i)| 2 ) < ^(|w(i')| 2 ) + -L||f || iao(Jifl) sup^(r)(t - t'), (17) 
for all t,t' G I, with t > t' . 

We also have the following a priori estimates for any solution w of the Navier-Stokes-a model in the sense 
of Definition 12.31 f adapted from gHl Corollary 3.2]): 

|w(i)| 2 < \w(t')\ 2 e-^-^ + J-^lfll 2 . (i - e-^-^), (18) 

Is 

1 /2 

Qf ||w(.s)|| W) < -^|w(f )| + A} /4 ^/ 2 (i - t') 1/2 , (19) 

1/2 

f T ll^w( S )|| 2 c(A) ,d S ) < ^i^|w(t')| 2 + ^M 2 + V V 2 \\' A M 2 (t - t>), (20) 

for all t', t E I, with £ > where c 2 is a universal constant and M 2 is a non-dimensional constant which depends 
only on non-dimensional combinations of the parameters v, \\ and |f ,h) ■ 

Again, we observe that if |w(f')| < R, for some R> Ro, where Ro is given by (jll[) . and for some t' G /, then 
it follows from (QJD that |w(t)| <R for all t>f. 

In Theorem 3.1], Vishik, Titi and Chcpyzhov proved the convergence of solutions of the Navier-Stokes- 
a model to solutions of the Navier-Stokes equations. This result was proved in the case where / G H and 
I = [0, oo), which was of interest to them. However, it is not hard to sec that the proof can be adapted to the 
case when / is an arbitrary interval and / G L°°(I,H). Since this result is going to play an important role in 
this article and for the reader's convenience we state it below: 

Theorem 2.5. Let L C K be an arbitrary interval and f G L°°(I,H). Let {w ra } be a bounded sequence in JFj 
such that each w„ is a solution of the Navier-Stokes-a n model on I , with a„ — > as n — > oo, and w„ — > w in 
t as n — y oo, for some w S Jj. Then w is a Leray-Hopf weak solution of the 3D Navier-Stokes equations on L. 
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2.5 Trajectory spaces 



The trajectory spaces made of solutions of the Navier-Stokes equations and of the Navier-Stokes-a model play an 
important role in the definition of the Vishik-Fursikov measure (see Definitions 12.41 and 12. 7|) since they connect 
this notion of solution with the corresponding equations. 

Let I C K be an arbitrary interval and R > 0. Consider the spaces Ci oc (I,H w ) and Ci oc (I, Bh(R) w ) 
endowed with the topology of uniform weak convergence on compact intervals in /. Then, Ci oc (I,H w ) is a 
separable Hausdorff locally convex topological vector space, hence completely regular, and C/ oc (7, Bh{R)w) is 
a Polish space, that is, a separable and completely metrizable space. 

For any interval J C /, we introduce the restriction operator 

n,/ : Cioc(I, H w ) — > Ci oc (J,H w ) 

u i ^ (nju)(t) = u(i), V* e J. 

It is clear that the restriction operator is continuous. Furthermore, if J is a closed subinterval of /, then IIj is 
also surjective and open. 

For each interval / in R and each t G I we define the projection operator U t by 

U t : Ci oc (I , H w ) — > H w 

u h-> II t u = u(t). 

which is also continuous, surjective and open. 

For the Navier-Stokes equations, we define the following trajectory spaces based on the Leray-Hopf weak 
solutions given by Definition 12. II 

Ui = {u G Ci oc (I; H w ) : u is a Leray-Hopf weak solution on /}, (21) 

U\ = {u <S Ci oc (I; H w ) : u is a Leray-Hopf weak solution on /}, (22) 

where / represents the interior of /. We endow these spaces with the topology inherited from Ci oc (I, H w ). 

The relation between the two spaces defined above is that u\ is the sequential closure of Ui with respect to 
the topology inherited from Ci oc (I; H w ). Furthermore it is clear that if I is open then Ui = u\. The difference 
appears when / is closed and bounded on the left, since we do not know, in general, whether the weak solutions 
in Jj\ are strongly continuous in H, from the right, at the left end point of the interval, as are those in Uj. 

Sometimes it will be useful to work with the following spaces 

Ui(R) = {u G Ci oc (I; Bh(R)w) : u is a Leray-Hopf weak solution on /}, (23) 

li\(R) = {u G Ci oc {I\ Bh(R) w ) : u is a Leray-Hopf weak solution on /}, (24) 

with the topology inherited from C; oc (J; H w ). 

As observed in if / is an interval closed and bounded on the left, then, for any sequence {Ri}°Z 1 of 
positive numbers with Ri > Rq, for all i G N, and Ri — > oo, we have the representation 

oo 

Uj = \Jul(Ri). (25) 

i=l 

And, if / is an interval open on the left, then, for any sequence {Ri]°^ 1 of positive numbers with Ri > Rq, for 
all and Ri — > oo and for any sequence { J n }^Li of compact intervals in / such that U^ =1 J n = I, we have 

the representation 
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u \= n u^xw (26) 

n— 1 i— 1 

As proved in [35], the spaces W/, Wj, Ui{R) and U\{R) are all Borel subsets of Ci oc (I, H w ) and, in particular, 
u\(R) is closed. 

For the Navier-Stokes-a model, we consider the solutions in the sense of Definition [531 Since these solutions 
belong to Ci oc (I; H), we define the following trajectory spaces 

Uf = {u £ Ci oc (I; H) : u is a solution of the Navicr-Stokcs-a model on /}, (27) 

Uf(R) = {u G Ci oc (I; Bh{R)) : u is a solution of the Navier-Stokes-a model on /}. (28) 

In order to compare with the solutions of the Navier-Stokes equations and since Ci oc (I,H) is included in 
Ci oc (I, H w ), we shall consider the spaces Uf and Uf{R) with the topology inherited from Ci oc (I, H w ). 

Here, we also have the same characterizations as the ones for the Navier-Stokes trajectory space. That is, 
if / is an interval closed and bounded on the left, then, for any sequence {Ri}°Zi °f positive numbers with 
Ri > Ro, for all i € N, and Ri — > oo, we have the representation 

oo 

U? = {JU?(Ri). (29) 

And, if / is an interval open on the left then, for any sequence {Ri}°Zi °f positive numbers with Ri > Rq, for 
all ieN, and Ri — > oo and for any sequence { J n }^ =1 of compact intervals in / such that U^ =1 J n = I, we have 
that 

oo oo 

"r = nU n l 1 W)- ( 3 °) 

The space Uf(R) is closed in the topology inherited from Ci oc (I, H w ) (and this implies that it is also closed in 
Ci oc (I, H) since H is continuously included in H w ). Indeed, since the solutions in Uf(R) are uniformly bounded 
by R in H, it suffices to show that Uf(R) is closed in Ci oc (I, Bh{R) w )- Since C[ oc (I, Bh(R) w ) is mctrizable, it 
suffices to work with sequences. Then, if {u„}„ is a sequence in Uf(R) which converges in Ci oc (I, Bh(R) w ) to 
an element u, then the a priori estimates ([T8)l . ([T9)l . ([20)l yield the compactness of this sequence in Ci oc (I,H w ) 
and in J 7 /. This compactness is sufficient to show that the limit function u is a solution of the Navicr-Stokcs-a 
model and, hence, belongs to Uf{R), proving that this space is closed. Since this space is closed, there is no 
need to consider spaces analogous to U\(R) and Uj, as done for the Navier-Stokes equations. 

Due to the representations ([29| and ([30)l . we see that Uf is an J>-set in Ci oc (I,H w ), in the case / is closed 
and bounded on the left, and it is an J-^-set, in the case I is open on the left. In any case, Uf is a Borel set. 

We now introduce an auxiliary functional space, yi, which is directly connected with the a priori estimates 
for the solutions of the Navier-Stokes equations and of the Navier-Stokes-a model, with suitable compactness 
property. First, let J be a compact interval in M, then we define 

y.j{R) = {xi£C loc {J;H w ) : \u(t)\ < R,\\u\\ L 2 {sAV) < -±^R + \\ fi vM(t - sf' 2 , and 

\\dM\ L ^ t , D{Ay) < - T7 f— R 2 + + u 5 / 2 xY 4 M(t -s),Va,tEj\, (31) 

A-j v ' A^ J 

where c = max{ci,C2} is a universal constant and M — max{Mi,M2} is a non-dimensional constant which 
depends only on non-dimensional combinations of the terms v : Ai and |jf .h) ■ With these choices of 
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constants, for any R > R , if \u \ < R, then the solutions of the Navier-Stokes equations and of the Navier- 
Stokes-a model with initial condition uq all satisfy the estimates in (|31[) . for subsequent times, which is possible 
thanks to the a priori estimates ©, © and (JTOJ) , and (jTSj , (HU), and (J20]). Thus, we have 

Uj(R), UJ{R) C yj(R), Vi? > R . (32) 

Now, let / be any interval in R and {J n }^Li be a sequence of compact intervals in / such that J n C J n +i 
and U^° =1 J„ = /. Consider also a sequence {-Ri}^ of increasing real number such that R\ > Rq and Ri — > oo 
as 2 — ?* cx). Define 

OO OO 

tt ; = n u n J>^(^)- ( 33 ) 

n— 1 ?—l 

Also, for a given R> R 0: we define 

oo 

y 7 (jz)= f|n^ Jn (j?). (34) 

n=l 

Observe that Uti !V/(-Ri) C J 7 / C Ci oc (I, H w ). The space yi(R) is independent of the choice of the intervals 

while the space yi is independent of the choice of both the intervals { J n }„ and the sequence 
although these properties are not really necessary. 

Lemma 2.2. Let J be a compact interval in K. Then, yj(R) is a compact subset of Ci oc (J, H w ). 

Proof. First, observe that since y.j(R) C C/ oc (J, Bh{R) w ) then y.j(R) is metrizable. Now, let {u„}„ be a 
sequence in y,j(R). It is clear that {u„} is bounded in L 2 (J, V) and {d t u n } is bounded in L 2 (J, D(A)'). 
Then, by Aubin's Compactness Theorem, we obtain that {u„} is relatively compact in L 2 ( J, H). Using all the 
information obtained before, we conclude that there exist a vector field u and a subsequence {u„ fc } such that 

u„ fc u in L°°(J;H); 
u„ fc ^ u in L 2 (J; V); 
d t u nk ^d t u mL 2 (J;D(AY); 
u nk -)• u in L 2 (J; H). 

Now, consider {w^}i a countable dense subset in D(A) (the existence of such a set follows from the fact that 
D(A) is separable). For each i G N, define the sequence {fl}k such that, for each k G N, fl(t) = (u nfe (t), w^) 
for all t g J. Note that {fl}k is a sequence of continuous functions from J to K, which is uniformly bounded 
and equicontinuous. Indeed, since 

< \t- s| 1/2 ||wJ D(A) ||<9 t u„J L 2 (M;ZW) , 

for all t,s£ J, then {/^}fc is equicontinuous. Thus, we can apply Arzela-Ascoli Theorem to obtain that {fl}k is 
relatively compact in C; oc ( J, R), for each ieN. By Cantor's diagonal argument, we can construct a subsequence 
{/£.} such that /|. converges to (u(-),w I ) as j —> oo, for all i 6 N. Now we use a triangulation argument to 
obtain the convergence of u„ fc . — >• u in C/ oc ( J, More precisely, given e > and v € H, there exist w G -D(A) 
such that |v — w| < ej (6R), and io G N such that ||w — w, \\d(A) < Aie/(6-R). By the convergence of {fl}j we 
conclude that there exists N G N such that sup tgJ |(u„ fc (t) — u(t),Wi)| < e/3, for all j > N and for all i G N. 



(^UnJ^jW^dT 
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Therefore, for all j > N, we have that 

sup|(u„ fe (t) -u(t),v) <sup|(u„ (t) -u(i),v-w)| 

tgj 3 tGJ 

+ sup|(u„ fe . (t) - u(t), w - w io )| + sup [(u„ fe . (i) - u(i),Wi )| 
teJ 3 teJ 3 

< SU P l u « fc , (*) ~ U WI(I V - w[ + |w - w 40 |) + | 
teJ 3 -J 

< 2i?|v - w| + 2R— ||w - w l0 ||_d (a) + £ < e. 

Ai z 

Finally, it is clear that u inherits the uniform estimates of the sequence in yj(R), so that u itself is in yj(R), 
completing the proof that yj(R) is compact. □ 



Lemma 2.3. Let I be any interval in R and R > 0. Then, the set yi(R) is compact in Ci oc (I,H w ). Moreover, 
yi(R) is metrizable. 

Proof. As observed in the end of Section 2.4 in [25] . compact subsets of Ci oc (I,H w ) can be characterized as 
the sets K for which, for every compact interval J C I, the subset HjK is equi-bounded with respect to the 
norm of H and equicontinuous with respect to the uniform structure of C; oc (J, H w ). Since Tljyi(R) C y,j(R) 
and yj(R) is compact as proved in Lemma 12.21 these conditions are met, and we have that yi(R) is compact. 
Furthermore, since yi(R) C Ci oc (I , Bh(R) w ) then yi{R) is metrizable. □ 



Lemma 2.4. Let I be any interval in R. Then, the space y^, endowed with the topology inherited from 
Ci oc (I,H w ), is a completely regular topological space. Moreover, yj contains the spaces U\ andUf. 

Proof. Since ^/ is a subspace of Ci oc (I, H w ) then 3^/ is a completely regular topological space. For the inclusions, 
suppose first that / is an interval closed and bounded on the left. Then, for every compact interval J C L 
containing the left end point of /, and for all R > Ro, where Rq is defined by pT|) . it follows from (pl2l that 
U\{R) C Tr, l yj{R). Therefore, using and (j2"5j) we conclude that u\ C yj. Now, in order to prove that 
Uf C yi, notice that from (|3"2"j) we have that Uf{R) C TYf l yi{R), for all R > Rq and any compact interval 
J C L containing the left end point of /. Thus, using (j34]l and (|29f we conclude that Uf C D^j- 

Now, suppose that / is open on the left. Again, we have from (p?2"j) that Uj{R) C yj(R), for all R > Rq and 
any compact interval J <Z I. Thus, by (|3"3")l and (f2l))) , we conclude that li\ C J/. Observe also that we have 
from ^ that Uj*(R) C yj{R), for all R > R and any compact interval J C I. Thus, ^ and ^ imply 
that Uf C yi . □ 



Remark 2.2. // / is an interval closed and bounded on the left and {Ri}i is a sequence of positive numbers 
with Ri > Rq, for all i G N, and Ri — > oo, then we have in fact showed in the proof of Lemma \2.J\ that 
Uf,Uj C Ui^i yi{Ri)- On the other hand, if L is an interval open on the left then Uf ,U\ C 3^/ but Uf and U\ 
might be not included in {J°^ 1 yi(Ri). 



Lemma 2.5. Let I be any interval in M. Then, the space yi(R) contains the spaces Uj(R) andUf(R), for all 
R > Rq, where Rq is defined by ill}) . 

Proof. It is clear from estimates ©, §9$ and (JlOj) that, if u e U$(R), for some R > Rq, then u e for 
all compact interval J C L. Thus, it follows from that U\{R) C yi(R). With an analogous argument we 
can also prove that Uf (R) C yi(R), for all R > R Q . □ 
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Next, we prove an important convergence result concerning the trajectory spaces U\(R) and Uf{R), based 
on Theorem 12.51 

Lemma 2.6. Let R > Rq, where Rq is defined by Ul)) . Let U\(R) and Uf(R) be given by \2J$ and i28\). 

respectively. Then, 

Kmdist yi(R) (U?(R),uK R )) = ° 

where 

disty l{H) {Uf{R),U\{R))= sup inf d(w,u), 

weWf(fl) ueujiR) 

and d is any compatible metric in yj (R) . 
Proof. Suppose by contradiction that 

\imdisty l(R) (Uf(R),Uj(R))^0. 
Thus, there exist e > and a sequence of {«„}„, with a n — > as n — > oo, such that 

sup inf d(w, u) > e, Vn £ N. 
yv£Uj n (R) uewf(R) 

Observe that, from the definition of the supremum, we have that, given r > 0, there exists w„ £ such 
that 

inf d(w n , u) > e — r, Vn £ N. 
uew|(fi) 

In particular, we can take r — e/2 and obtain 

£ 

inf d(w„,u) > -, Vn £ N, 

so that 

rf(w n ,u) > |, Vn £N and Vu eU\{R). (35) 

On the other hand, we have that |w n (£)|# < R for all t £ I and for all n G N. Thus, estimates (fT9|) and 
([2^| imply that {w^jn is bounded in J| (see ©). Therefore, there exist a subsequence {w ni }j of {w n }„ and 
a function u £ J-\ such that w, i; — > u with respect to r. Using Theorem 12. 5[ we conclude that u £ U\{R). 
Moreover, for any compact interval J C /, we have that {w„}„ is in yj(R). Then using Lemma l2.2l we conclude 
that there exists a subsequence of {w„}„ that converges in the topology of weak converge in H uniformly on 
J. We can now use Cantor's diagonal argument to obtain a subsequence {w n) , }„ fc that converges to u in the 
topology of weak converge in H uniformly on J, for any compact interval J <Z I. Observe that this lead us to 
a contradiction with (|35|) since d is a compatible metric with the topology of yi(R). □ 

2.6 Statistical solutions 

The notion of statistical solutions that is considered here was introduced by Foias, Rosa and Temam, in [531 HI] 
(see also [3SJ US]). We recall this definition in the context of the Navier-Stokes equations and introduce a 
corresponding definition for the Navicr-Stokcs-a model. 
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2.6.1 Time-dependent statistical solutions 

We start with the definition of Vishik-Fursikov measure for the Navier-Stokes equations which will give rise to 
the statistical solutions for the Navier-Stokes equations. 

Definition 2.4. Let I C R be an interval. We say that a Borel probability measure p in Ci oc {I,H w ) is a 
Vishik-Fursikov measure over I if p satisfies the following 

(i) p is carried byU\; 

(ii) tH>/ M ,|u(t)| a dp(u) €L% C (I); 

(Hi) if I is closed and bounded on the left, with left end point to, then for all tp £ \t we have that 

lim / ri(|u(f)| 2 )dp(u) = [ M\n{U)\ 2 )dp{n), 

where * := {tp £ C 1 ([0, oo)) : tp > 0, tp' > and sup t>0 tp' '(t) < oo}. 

Observe that in this definition we only require the measure p to be carried by li\ . But we really want to have 
p carried by Ui, which is the trajectory space of Leray-Hopf weak solutions. And this is what in fact happens. 
More precisely, in Theorem 4.1], it was proved that for an arbitrary interval I Cl, any Vishik-Fursikov 
measure over I is carried by U[. 

Next we present the definition of a Foias-Prodi statistical solution, which is a family of measures on the 
phase space, satisfying a Liouville-type equation and some regularity properties. Let us denote by F the operator 
defined on V as F(u) = f — vAvl — B(u, u), with values in V' . A function $ : H — > R is called a cylindric test 
function if 

$(u) = ^((u,vi),...,(u,v fc )), 

where k G N, ip is a continuously diffcrcntiable real- valued function on R fe with compact support, and Vi, . . . , 
belong to V. For such $, we denote by its Frechet derivative in H, which has the form 

fc 

= Yl 9 M( u ^i )>•••> (u. v fe))vj, 

3=1 

where djtp is the derivative of tp with respect to its j-ih coordinate. 

Definition 2.5. A family {pt}t>o of Borel probabilities on H is a statistical solution of the 3D Navier-Stokes 
equations if it satisfies: 

(i) the Liouville type equation 

j t J $(u) dp t (u) = J <F(u), $'(u)) du t (u), 

in the distributional sense in t > 0, for all cylindric test functions $; 
(ii) the function 

1 1 y / ^(u) dnt(u) 
Jh 

is measurable in t > for all continuous functional <p '■ H — > R; 
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(Hi) the function 

t i-> / |u|^- dfH(u) 

JH 

belongs to Lj^ c (0,oo); 
(iv) the function 

tH. f \\u\$dnt(yL) 

JH 

belongs to Lj oc (0,oo); 
(v) the mean strengthened energy ineguality holds, i.e., 

^ dfitiu) +v Sj' (|u| ^ ||u||2y dA<t(u) - j H ^ (|u| ^ (f ' u) ^ t(u) 

in the distributional sense in t > 0, for all ip G 
(vi) and the function 

1 1 y I ip(\u\ 2 H ) d/it(u) 

JH 

is continuous at t = 0, for all ip G \l/ . 

In [25] . it was proved that for any Borel probability measure uq on H such that J H \u\ 2 dpo(u) < oo, there 
exists a Vishik-Fursikov measure p over I = [to, oo) such that Tlt p = fio- Furthermore, if p is a Vishik-Fursikov 
measure then fit := ILtPi t G /, is a statistical solution in the sense of Definition 12.51 This yields a particular 
type of statistical solution: 

Definition 2.6. Let I C be an arbitrary interval. A Vishik-Fursikov statistical solution of the Navier- 
Stokes equations over I is a statistical solution {pt}tei such that p t = HtP, for all t G I, for some Vishik-Fursikov 
measure p over the interval I. 

Inspired by the definition of a Vishik-Fursikov measure, we define the a- Vishik-Fursikov measure which will 
give rise to the statistical solutions for the Navier-Stokes-a model. 

Definition 2.7. Let L C K be an interval and a > 0. We say that a Borel probability measure p a in Ci oc {L, H) 
is an a-Vishik-Fursikov measure over L if p a satisfies the following 

(i) p a is carried byUf; 

(ii) t^f u? \w(t)\ 2 d Pa (w)eL% c (I). 

Since wc chose to work with solutions of the Navier-Stokes-a model in the sense of Definition [531 an d thanks 
to condition (iii) of that definition, it is natural to define a- Vishik-Fursikov measures on C; oc (J, H). Since lAf 
is contained in Ci oc (L, H) as a set, the condition (ii) of the Definition 12.71 above makes sense. 

As another remark, since C; oc (J, H) is continuously included in C/ oc (/, H w ) and the space lAf(R) is closed in 
Ci oc (L, H w ) (see Section [5T5]) . there is no need to consider a space analogous to u\, as in Definition 12.41 

Remark 2.3. Observe that Ci oc (L,H) is a Polish space (since it is Frechet and separable) so that every Borel 
probability measure p is tight. Moreover, since the Borel sets of H are the same as the Borel sets of H w ( see 
\2iy ) then the Borel sets ofCi oc (I,H) are the same as the Borel sets of Ci oc (I,H w ). Therefore, every Borel 
probability measure in Ci oc (I , H w ) is tight. 
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We also define an a- Vishik-Fursikov statistical solution of the Navier-Stokes-a model over an arbitrary 
interval / as a family {pf} t< zi of Borel probability measures on H such that pf = Htp a , for all t G /, for some 
a- Vishik-Fursikov measure p a over the interval I. 

The existence of a- Vishik-Fursikov measures is easy to obtain. For instance, any Dirac measure in the 
trajectory space Uf is an a-Vishik-Fursikov measure over i\ Furthermore, given an initial Borel probability 
measure po in H with finite energy we can construct an a- Vishik-Fursikov measure p a over [0, oo) such that 
IloPa = Mo- Indeed, since the Navier- Stokes- a model is well-posed, the solution semigroup {S a (t)}t>o is well- 
defined. Moreover, the operator S a (-) : H — > Ci oc (I,H), defined as S a (-)wo = w(-) for wo G H, where 
S a (t)\VQ = w(£) for all t G /, is continuous. Therefore, given an initial Borel probability measure po on H, we 
may define the Borel probability measure p a as 

p a {E) = poiSai-^E), for all E G C loc (I,H) Borel, where I = [0,oo). 

By construction, it is clear that p a is carried by the set Uf. Moreover, since for each w G Uf it holds that 
|w(i)| < |w(0)| + 1/(A 2 (■)po, then, using the Change of Variables 

Theorem (see Section [22]), we obtain that 

/ |w(t)| 2 dp Q (w) < / |w(0)| 2 dp Q (w) + -l T ||f||| oo(/!ff) 
JUf JUf A l v l 

= [ |(5 Q (.)u)(0)| 2 d M o(u) + -^||f||i. (/ , H) 
= I |u| 2 ^ (u) + -^||f||^ (7iff) . 

JH A \ v \ 

Consequently, since po has finite energy we obtain that t h- » J ua \w(t)\ 2 dp a (w) G Lf£ c ([0, oo)). Therefore, p a is 
an a-Vishik-Fursikov measure, and it is straightforward that no/?a — Mo- 

Nevertheless, since the well-posedncss for the 3D Navier-Stokcs equations is an open problem, the abstract 
definition of Vishik-Fursikov measure is essential in the context of this article. 

2.6.2 Stationary statistical solutions 

As already mentioned in the Introduction, stationary statistical solutions are valuable in the study of turbulence 
in statistical equilibrium in time, yielding, in particular, rigorous proofs of important statistical estimates. The 
concept of stationary statistical solution represents a generalization of invariant measures for the semigroup 
generated by an equation. For instance, in 2D, since the Navier-Stokes equations has a well-defined semigroup, 
stationary statistical solutions are, under certain hypothesis, precisely the invariant measures for the semigroup; 
see [3T] for more details. 

Suppose / is an interval unbounded on the right, hence having one of the following forms: I = [to,oo), 
/ = (to, oo) or I = R. We introduce the time-shift operator a T defined for any r > by 

&t '■ Ci oc (I,H w ) —> Ci oc (I,H w ) 

U I y (<T T U)(t) = U(t + T), W G i". 

An invariant Vishik-Fursikov measure over / is a Vishik-Fursikov measure p which is invariant with respect 
to the translation semigroup {cr T } r >o, in the sense that a T p = p for all r > 0, i.e., p(E) = p(a~ 1 E), for all 
Borel set E in Ci oc (I,H w ). 

The family of projections {H t p} t £i of an invariant Vishik-Fursikov measure p has the property that any 
statistical information 

/ #u)dCM(u) = / 0(u(f))dp(u) 

JH JH 
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is independent of the time variable t £ I, for any (f> £ Cb{H w ). In fact, the measure ILp itself is independent of 
t and is a statistical solution in the sense of Foias-Prodi which is time-independent (what is called a stationary 
statistical solution in the sense of Foias-Prodi). This yields a particular type of stationary statistical solution 
as stated precisely below. 

Definition 2.8. A stationary Vishik-Fursikov statistical solution on H is a Borel probability measure 
po on H which is a projection po — TltP, at an arbitrary time t £ I , of an invariant Vishik-Fursikov measure p 
over an interval I unbounded on the right. 

3 Convergence of statistical solutions of the a-model as a vanishes 

In this section we present the main results of this paper. We prove that, under certain conditions, statistical 
solutions of the Navier-Stokes-a model converge to statistical solutions of the 3D Navier-Stokes equations as 
ol — y 0. We first prove this result for time-dependent statistical solutions and then we address the particular 
case of stationary statistical solutions. 

As defined in Section 12721 we mean by p a — 1 p in V(X) that 



Observe that, by Theorem 12.21 in the case when X is a completely regular Hausdorff space, which is the case 
for us, this is exactly the weak convergence — » discussed in Section l2~2l 

3.1 Time-dependent statistical solutions 

In this section we state some results concerning the convergence of a- Vishik-Fursikov measures and statistical 
solutions as a vanishes. We will see that it suffices to impose some condition of uniform boundedness on the 
mean kinetic energy over the a- Vishik-Fursikov measures (or, equivalently, on the a- Vishik-Fursikov statistical 
solutions). This will assure the tightness of the a- Vishik-Fursikov measures and yield a convergent subsequence. 
The convergence of these measures will imply the convergence of the associated statistical solutions in a suitable 
sense, described as follows. 

Let J C M be an arbitrary interval. We say that a family {{pf}tei}a>o of a- Vishik-Fursikov statistical 
solutions over I converges, as a — » 0, to a Vishik-Fursikov statistical solution over / if there exists a Vishik- 
Fursikov statistical solution {p, t }tei over I such that 



Recall that a Vishik-Fursikov statistical solution (a- Vishik-Fursikov statistical solution) over an interval I is 
a family of Borel probability measures {p t }t£i given by pt = n t/ o, for all t £ 1 . where p is some Vishik-Fursikov 
measure (a- Vishik-Fursikov measure) over I. 

Thus, if we have a family of a- Vishik-Fursikov measures {p Q } Q >o that converges to a Vishik-Fursikov measure 



then, as a simple consequence of the Change of Variables Theorem (see Section |2"?2"]) . we obtain the convergence 
of the a- Vishik-Fursikov statistical solutions {Htp a }t£i to the Vishik-Fursikov statistical solution {n t p} tg j. 

The measures p and p a are defined on Ci oc (I , H w ) and Ci oc (I , H), respectively. Since Ci oc (I,H) is included 
in Ci oc (I, H w ) we may consider them as measures on Ci oc (I, H w ). In fact, since they are actually carried by the 





(36) 



P in P(yi), i.e., 




(37) 
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space of solutions Ui and Uf , respectively, and these spaces are included in yj, we may consider them restricted 
to yj. With this in mind and for the sake of simplicity we consider, in what follows, the measures p and p a as 
measures on yj. This is the same reason why we used in the convergence f|3T[> . 

We state first a lemma which is essentially Theorem 12.21 translated into the framework of interest in this 
section. 

Lemma 3.1. Let I be an arbitrary interval in R and {p a }a>o be a family in V(yi,t). Suppose that for 
every e > there exists R > such that p a {yi \ yi(R)) < £• for all a > 0. Then, there exists a sequence 
{Pa„}n C {p a }a>o, with a n — > as n — > oo , such that p an p in V(yi) as n — > oo. 

Proof. First, observe that since {a G R : a > 0} is a totally ordered set we can subtract a sequence {p an }n from 
{Pa}a, with a n — > as n — > oo. Moreover, observe that yj is a completely regular Hausdorff space and yi{R) 
is compact, for any R > 0. Therefore, it is clear that the sequence {p an }n<£N fulfills the hypothesis of Theorem 
12.21 Therefore, there exist p G V{yi,t) and a subsequence, which we still denote by {p an }n C {p a }a>o, such 
that p an A p in V{yi). □ 

We would like to highlight that the use of a sequence in Lemma [57X1 is only for the sake of simplicity. Notice 
that Theorem 1 2. 2 1 can be applied to the net {p a } a , as in Lemma l3~Tt and it yields the existence of a convergent 
subnet of {p a } a - Then, all the subsequent results are also valid if we work with a convergent subnet instead of 
a sequence. 

Observe that in order to apply the last lemma to a family of a-Vishik Fursikov measures we need the 
sequence to be uniformly tight. A natural condition that yields this uniform tightness is obtained by imposing 
a uniform boundedness condition on the mean kinetic energy of this family. This uniform boundedness of the 
mean kinetic energy will also be important to yield that the limit measure has finite mean kinetic energy and 
is a Vishik-Fursikov measure. This uniform boundedness can be imposed in different ways. We start with the 
following: 

Proposition 3.1. Let L be any interval in R. Let {p a }a>o be a family of Borel probability measures on yi 
such that, for each a > 0, p a is an a-Vishik-Fursikov measure over L and sup tg j Jy |w(t)| 2 dp Q (w) < C , for all 
a > 0, for some constant C > 0. Then, there exists a sequence {p a „}n C {p a }a>o, with a n — > 0, converging to 
a Borel probability measure p in yj . Moreover, p is carried by U\ and HjP is a Vishik-Fursikov measure over L. 

Proof. In order to prove the convergence we will use Lemma 13.11 First, recall that every a- Vishik-Fursikov 
measure is tight ( Remark 12. 3[) . Next, we check the uniform tightness condition. Observe that it is enough to 
prove that for all e > there exists R > Rq such that p a (Uf \ Uf(R)) < e for all a > 0. Indeed, since p a is 
carried by Uf we have that 



where the inequality holds since Uf(R) C yj(R), which implies that yi(R) c C {Uf(R)) c . 

Let R > R a . If Uf \ Uf{R) is not empty, then for all w G Uf \ Uf(R), there exists t w G / such that 
|w(t w )| 2 > R. Thus, sup te/ |w(i)| 2 > R and 



p a {yi \ yi(R)) = p a {yi n (yi(R)) c ) = p a (uf n (yi(R)) c ) 

< Pa (Uf n (Uf(R)) c ) = Pa (Uf\Uf(R)) : 




(38) 



Since for all w G Uf we have that 



|w(i)| 2 < |w(t')| 2 e 



Aii/(t-t') 



1 



fll? 




Xii/(t-t') 



+ 



\\v 2 



L°°(t',t;H) 
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for all t,t' £ I with t > t', then 

[ sup|w(i)| 2 dp Q < / \™(t')\ 2 dp a + -^\\{\\ 2 Lao{IH) . 

JUf t>t> Juf A l v 

Therefore, by the last inequality and the hypothesis sup tg/ Jy \w(t)\ 2 dp a (w) < C, for all a > 0, it follows that 

/ sup \w(t)\ 2 dp a < Ci, 
JUf t>t< 

for all t' el and for all a > 0, where C x = C + l/(A?i/ 2 )||f |||oc (/iiJ) - 

Also observe that, for all w G Uf, w G Ci oc (I, H) so that the function defined by f(t') = sup t>t / |w(i)| 2 for 
t' G J, is continuous. Thus, for any to € I and any sequence {t' k }k C / such that t' k — > t^, using the Monotone 
Convergence Theorem, we find that 

/ sup \w(t)\ 2 dp a < / \w(t')\ 2 dp a < Ci, Vi G I and Va > 0. 

JUf t>t Juf 



In particular, we obtain 

/ sup \w(t)\ 2 dp a < d, Ma > 0. 
Juf tei 

We use the last estimate in ([55)1 to obtain, in the case that Uf \Uf(R) is not empty, that 

p a {Uf\Uf{R))<^, Va>0. 

If Uf \ Uf(R) is empty, then this estimate is trivially valid. Then, given e > 0, take R = max{2Ci/e, Rq}, so 
that 

Pa (Uf\Uf{R))<£, Va>0. 

This shows that {p a } Q >o is uniformly tight. 

Then, we can apply Lemma l3~Tl and obtain the existence of a measure p G Vtyi, t) and a sequence {p Q „}n C 
{/3q}q>0i with a ra — )■ as n — > oo, such that p Qn — ^ p in Vtyi) as n — s> oo. 

Next, we prove that p is carried by u\. In order to do so, we define, for each e > and i? > i?o, the set 

X(i?) = {u G : &sty l{R) {u,U\{R)) < £}• 



Observe that yi(R) \y e (R) and 3^/2(-R) are disjoint closed sets in yi(R). Since yi(R) is a compact Hausdorff 
space then, by Urysohn's Lemma (see e.g. [T]), for each e > 0, there exists a continuous function ^ : yi(R) —> 
[0,1], such that <pf(u) = 1 for all u G y s/2 (R) and pf(u) = for all u G ^j(iJ) \ ^ e (i?). Now, we define 
an extension of ipf to IVj, ^ : — > [0, 1] as (pf (u) = ipf-(u), for all u G yi{R) and </5f (u) = 0, for all 
u G yi \ yi(R). Since yi{R) is closed, it is easy to see that (pf is upper semicontinuous. Then, using that 
Pa n ~ ^ P m "P{yi) and Lemma 12. 11 we obtain that 



/ tpf (u)dp(u) > limsup / <pf(u)d Pan (u) = limsup / ^(u)d/o an (u). 



(39) 



Since we have, by Lemma l2~6l that limn-^oo dist y T rfy (Uf™ (R) ,U\(R)) = then, given e > 0, there exists iVo G N 
such that, for each m > Nq, 

disty l(R) (Uf™^(R),uj(R)) <e/2, 
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which implies that for all u G U" m (R), m > N , diaty I t R \(u,Uj(R)) < e/2. In other words, given e > there 
exists N Q e N such that Uj m (R) C 3> e / 2 (-R), for all m > N Q and therefore <^f (u) = 1 for all u G U? m (R), for 
all m>N . Thus, for all m>N , 



f <p?(u)d P a ™ (u) = p Qm (Z#» (i?)) > 1 - 
hence, using the last estimate in we conclude that 

J $i?(u)dp(u)>l-^. 



Therefore, 



p(y e (ie)) > / ^f(u)dp(u) > i - 



Since U\(R) = n"^^ (-R) for an y sequence of positive numbers ej — > 0, we obtain that 

> 1-d/R, Vi?>i? . 

Now, in order to prove that p{U\) = 1, we first suppose that / is closed and bounded on the left. In this case we 
can write U\ = yj°^L l bl\{Ri) , for any sequence {Ri}i with Rq < Ri < R4+1, for all i G N, and i?i —5- 00. Observe 
that u\{Ri) C U\(R i+ i), for all i G N. Therefore, 

P {U\) = p (\JUi(Ri^j = Um > Um (l " ^) = 

Otherwise, we can write Wj = H^Li Ui=i nj^Wj C^*)> f° r an y sequence with Ro < Ri < Ri+i, for 

all i £ N, ft -K» and for any sequence {J„}„ of compact subintervals of / such that J„ C J n +i, for all 
?i G N, and U„J„ = /. Observe that, for all compact interval J C / it is true that U\(R) C IIJ 1 t/j(i?) so that 
p(IIj 1 Wj(i?)) > 1 — C\/R. Moreover, it is easy to see the monotonicity properties 

00 00 

U n^X U n ^l +1 

i=l i=l 

for all n G N, and 

irj^ c ) , 

for all ii,ieN. Therefore, 



(00 00 \ / 00 \ 

n u u j1 u i = p u n ^ iw i (^) 
n=\i=\ ) " °° \i=l / 

= lim lim ^(nT 1 ^* (iJi)) > lim lim (l - — ) = 1. 

Thus, p(w|) = 1. 

It remains to prove that i i-)- J!y |u(7j)| 2 <i/9(u) belongs to Lf£ c (I). In that direction, we define an increasing 
sequence of cut-off functions {4>m}m, that is, for all M > we have that G C°°([0, 00)), < 4>m < 1, 
<?W(f) = 1 for < r < M, </>(r) = for r > 2M and <Pm < </>a/+i- Now, take any t' £ I and observe that the 
function fM,k, defined by /M,fc(u) = 4>m (\PkU-(t')\ 2 )\PkU.(t')\ 2 , where Pk is the Galcrkin projector (see Section 
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12. lp , belongs to Cf>(3^) and fM,k( u ) — I u (*')| 2 j f° r au u e Then, from the convergence of p an to p together 
with the hypothesis sup tg/ fy \w(t)\ 2 dp an (w) < C, for all n € N, we obtain 

/ fM,k(vL)dp(u) = lim / / M ,fc(u)dp Qn (u) < limsup / |u(i')| 2 <i/3a„ (u) < C. 

Jyj "^^Jyj n^co Jy T 

We can pass to the limit in the last inequality as M —> oo and, using the Monotone Convergence Theorem, 
we obtain that 

/ \P k u(t')\ 2 dp(u) = lim / f M , k (u)dp{u) < C. 

hi M ^°° hi 

Again, using the Monotone Convergence Theorem, we can pass to the limit as k — > oo to find that 

/ \u(t')\ 2 dp{u) = lim / \P k u(t')\ 2 dp{u) < C. 

hi k ^°° hi 

Now, since t' S I is arbitrary we obtain that t fy \u(t)\ 2 dp(u) belongs to L°°(I). □ 
The previous result has a corresponding statement in terms of statistical solutions, which we write as follows. 



Proposition 3.2. Let I be any interval in R and let {{p-f}tei}a>o be a family of a-Vishik-Fursikov statistical 
solutions over I, such that sup tg/ f H |w| 2 d^"(w) < C , for all a > 0, for some constant C > 0. Then, there 
exists a sequence {{fJ^ n }tei}n C {{p't}t<=i}a>o, with a n — > 0, converging to a Vishik-Fursikov statistical solution 

Proof. By definition, for each a > 0, there exists an a- Vishik-Fursikov measure p a over / such that pf = Htp a , 
for all t £ I. Observe that the family {p a } a >o fulfills the hypothesis of Proposition 13.11 then there exists a 
sequence {p a „}n C {p a }a>o, with a n — > 0, that converges to a Borcl probability measure p in yj such that 
HjP is a Vishik-Fursikov measure over I. Thus, as a simple consequence of the Change of Variables Theorem 
(see Section |2~2]) . we obtain the convergence of the a n -Vishik-Fursikov statistical solutions {Ii t p an }tei to the 
Vishik-Fursikov statistical solution {H t p} te f. In other words, {{Pt n }tei}k converges to the Vishik-Fursikov 
statistical solution {fit} te f, where pt = n t p, for all t € I. □ 

A different way of obtaining the uniform tightness condition in the family of a- Vishik-Fursikov measure is 
to assume that the interval / is bounded and closed on the left and impose uniform boundedncss on the initial 
mean kinetic energy over the a- Vishik-Fursikov measures. This is done in the next result: 

Corollary 3.1. Let L be an interval in R which is bounded and closed on the left, with left end point to. Let 
{pa}a>o be a family of a- Vishik-Fursikov measures over I such that fy \w(to)\ 2 dp a (\v) < C , for all a > 0, for 
some constant C > 0. Then, there exists a sequence {pa n } n C {p a }a>0; with a n — > 0, that converges to a Borel 
probability measure p in yj such that p is carried by U\ and II^p is a Vishik-Fursikov measure over L. 

Proof. We only need to check that the family {p a } fulfills the hypothesis of Proposition l3.il In order to do so, 
observe that for all w S Uf and for all i G I with t > to, 

|w(t)| 2 < |w(t )| 2 e-^^) + -IjHf \\U {t(> _ t;ff) (l - e-^'-'o)), 



therefore 



/ \w(t)\ 2 d Pa (w) < [ \Mto)\ 2 dp a (w) + -^\\t\\U {to , t . H) 

JUf JUf 
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Since by hypothesis Jy \w(to)\ 2 dp a (w) < C, for all a > 0, and p a is carried by Uf, then 

sup sup / \w(t)\ 2 dp a (u) < Ci, 
q>o tei Jy! 

where d = C + I/^iWl^hy □ 

As before, the previous result has a corresponding statement in terms of statistical solutions, which wc write 
as follows. 

Corollary 3.2. Let I be any interval in M. and let {{/i"} te /} Q ,>o be a family of a-Vishik-Fursikov statistical 
solutions over I, such that J„ |w| 2 d/i" o (w) < C, for all a > 0, for some constant C > 0. Then, there exists 
a sequence {{Pt"}tei}n C {{pf}tiEi}a>o, with a n — > 0, that converges to a Vishik-Fursikov statistical solution 

In the last two results we obtain a limit of a- Vishik-Fursikov measures over an interval / which is a Vishik- 
Fursikov measure only in the interior of the interval I. If / is an interval open on the left then the limit is in 
fact a Vishik-Fursikov measure over the whole interval /. The problem is when / is closed and bounded on the 
left, since in this case we may lose the continuity at the left end point of / of the strengthened mean kinetic 
energy for the limit measure (see condition (Hi) of the Definition 12. 4[) . This problem is, in fact, inherited from 
an analogous problem for individual weak solutions, as described in Section |2~51 and which led us to introduce 
the spaces u\ (R) and u\ . 

In the next result we impose some conditions over a family of a- Vishik-Fursikov measures in order to have 
that the limit is in fact a Vishik-Fursikov measure over the whole interval /. 

In fact we impose conditions only on the initial measures II to p a , where to is the left end point of the interval 
/. These initial measures shall converge in a slightly stronger sense to a measure fiQ in H w . Since the er-algebra 
of Borel sets with respect to the weak topology of H coincides with that with respect to the strong topology 
we may consider fj,Q as a Borel probability either on H or on H w . 

Theorem 3.1. Let L be an interval in R bounded and closed on the left, with left end point to, and let {p a }a>o 
be a family of a-Vishik-Fursikov measures over L. In addition, suppose that there exists a Borel probability 
measure fiQ on FL such that 

ft) Ih l u | 2 <Wu) < oo; 
(ii) Tito Pa P-o in V{H W ); 
(Hi) and for all tp G \J r 

lira / il>{\u(t )\ 2 )d Pa (u) = [ ^(|u| 2 )^ (u). 

Then, there exists a sequence {p an }n C {pa} Q >o, with a n — >• 0, that converges to a Borel probability measure p 
in yj such that p is a Vishik-Fursikov measure over I and Tlt P = Pa- 
Proof. Observe that by taking ijj = 1 we can sec that the family {p a } a >o fulfills the hypothesis of Corollary [XT] 
Thus, there exist a sequence {p a „}n C {p a }a>o, with a„ — > 0, and Borel probability measure p in yj such that 
Pa n — 1 P in V{yi), as n —¥ oo, p is carried by U '|, and Hjp is a Vishik-Fursikov measure over I. Then, in order 
to obtain that p is a Vishik-Fursikov measure over /, it remains to prove that for all ip S ^ it holds true that 

hm / ^(|u(t)| 2 )dp(u) = / V(|u(t )| 2 )dp(u). (40) 
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Observe that for any <f> G Cb{H w ) the function defined by ip = <f) ° Iltota belongs to <p G Cb(yi). Then, it is 
straightforward from the convergence of p an , from hypothesis (ii) and from the Change of Variables Theorem 
that 

/ 0(u)djti O (u) = / 0(u(t o ))d/o(u), for all fj> G C b (H w ). 
jh Jyi 

Therefore, since if™ is a completely regular Hausdorff space and /io,n to p € V(H w ;t) we obtain that n to p = fio 
(see ©). 

Now, take ip G ^ and observe that, since any u G is weakly continuous at to and ^ is nondecreasing and 
continuous, then -0(|u(io)| 2 ) < lim inf t->t VKi 11 ^)! 2 )- Therefore, it follows from Fatou's Lemma that 

if)(\u(t )\ 2 )dp(u) <liminf f ifj(\u(t)\ 2 )dp(u). (41) 

Let </)Af be a function in C 1 ([0,oo)) defined as 4>m{t) = 1, for all < r < M, <^>m(»") = 0, for all r > 2M, 
<j>M < 0m+i, and < 0jv/ < 1, for all M G N. It is clear that ip m (u) := ip(\P m u(t)\ 2 )(f>M (\Pmu(t)\ 2 ) belongs to 
Cb(yi), for any t G /, so that 



^(|P m u(t)H0 M (|P m u(i)| 2 )dp(u) = lim / VXI^™uWr)^(l-PmuWr)rfpa„(u)- 

Since < 1 and if>(\P m u(t)\ 2 ) < if>(\u(t)\ 2 ) we find that the right hand side of the last identity is bounded 
above by lim sup^^^ fy if)(\u(t)\ 2 )dp ari (u). On the other hand, using the Monotone Convergence Theorem 
twice we obtain that 

lim lim / ^{\P m VL{t)\ 2 )M\P m n{t)\ 2 )dp(u) = [ 4>(\u(t)\ 2 )dp(u). 

m— yoo M— - - 



Therefore, for any t G /, 



yi Jyi 



if)(\u(t)\ z )dp(u) < limsup / iP(\u(t)\ 2 )dp an (u). (42) 

Vi n->oo Jyj 



Now, using inequality (|32"|). estimate (TIT]) and the facts that p a is carried by Uf and ip is nondecreasing, we 
obtain 

limsup / tf)(\u(t)\ 2 )dp(u) < lim sup lim sup / ip(\u(t)\ 2 )dp an (u) 
t—>to Jy T t— H n— >oo jy 7 

< limsup limsup ( / V(|u(t )| 2 )dp Qn (u) + - — ||/||£°o(i,m supip'(r)(t - to)) 

t->t n^oo VJ^j Alt' r >0 / 

< limsup / V(|u(£o)| 2 )d/9 Q „(u) < / V(|u| 2 )dMo(u), 

where the last inequality follows from hypothesis (Hi). 

To conclude, we use that H to p = po in the last inequality together with (|4"Tj) to obtain that (|4U)) holds 
true. □ 

Now, we write the previous result in terms of statistical solutions. 

Theorem 3.2. Let I be an interval in M. bounded and closed on the left, with left end point to, and let 
{{p't}tei}a>o be a family of a-Vishik-Fursikov statistical solutions over I. In addition, suppose that there 
exists a Borel probability measure po on H such that 

(i) Ih l u l 2 <Wu) < oo; 
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(H) Mt* A Mo in V{H W ); 
(Hi) and for all ip G ^ 

lim / V>(|u| 2 >K(u) = f ^(|u| 2 )d/i (u). 

Then, there exists a sequence {{/i""} te /}„ C {{pf}t<£i} a >o> a„ — > 0, that converges to a Vishik-Fursikov 
statistical solution {p t }tei such that p ta = Po- 

Proof. It follows from the definition of a- Vishik-Fursikov statistical solution the existence of an a- Vishik- 
Fursikov measure p a over / such that p" = lit Pa, for all t € I. It is clear that the family {p a } a >o fulfills 
the hypothesis of Theorem 13.11 so that there exists a sequence {p a „}n C {p a }a>o that converges to a Vishik- 
Fursikov measure p over /. Therefore, the corresponding sequence {{Pt n }tei}n converges to the Vishik-Fursikov 
statistical solution {p t }t£i 7 where pt = II t p, for all t £ I . □ 

A case of particular interest for approximation purposes is when the measures p a have the same initial 
projection Yl to p a = po, for all a. This leads us to the following two corollaries, in terms of Vishik-Fursikov 
measures and Vishik-Fursikov statistical solutions, respectively. 

Corollary 3.3. Let I be an interval in R closed and bounded on the left, with left end point to, and po a Borel 
probability on H such that j H \u\ 2 dpo(u) < oo. Let {p a }a>o be a family of a-Vishik-Fursikov measures over 
I such that Ht p a = Po, for all a > 0. Then, there exists a sequence {p a „}n C {pa}a>cb with a n — > 0, that 
converges to a Borel probability measure p in yj such that p is a Vishik-Fursikov measure over L and II toJ o = p,Q. 



Corollary 3.4. Let I be an interval in M. closed and bounded on the left, with left end point to, and p,Q a Borel 
probability on H such that J H \u\ 2 dfio(u) < oo. Let {{^"}te/} Q >o be a family of a-Vishik-Fursikov statistical 
solutions over L such that pf = po, for all a > 0. Then, there exists a sequence {{H t p an }tei}n C {{p?}tei}a>o, 
with a n — > 0, that converges to a Vishik-Fursikov statistical solution {p t }tei such that pt = po- 

3.2 Stationary statistical solution 

We say that a family {p a } a >o of stationary a- Vishik-Fursikov statistical solutions converges to a stationary 
Vishik-Fursikov statistical solution if there exists a Borel probability measure p on H such that p is a stationary 
Vishik-Fursikov statistical solution and 

lim / cp(u)dp a (u) = / 0(u)d/i(u), 

for all (f> G C b (H w ). 

Theorem 3.3. Let {p a } a >o be a family such that, for each a > 0, p a is a stationary a- Vishik-Fursikov 
statistical solution. Suppose that J H |u| 2 d^ Q (u) < C, for all a > 0, for some C > 0. Then, there exists a 
sequence {p a „} n C {p a } a >o, with a n —> 0, that converges to a Borel probability measure p on H. Moreover, p 
is a stationary Vishik-Fursikov statistical solution. 

Proof. Since p a is a stationary a- Vishik-Fursikov statistical solution, for each a > 0, there exists an invariant 
Vishik-Fursikov measure p a over an interval / unbounded on the right, such that p a = n t p Q , at any time tel. 
Hence, using the Change of Variables Theorem it follows that 



f \u{t)\ 2 dp a (u) = [ |u| 2 d^(u), Vt G I, Va > 0. 
Jyj J h 



2G 



From the hypothesis j„ \u\ 2 dp a (u) < C, for all a > 0, and the last identity, we obtain that 

sup / \u{t)\ 2 dp a (u) < C, Va > 0. 

Thus, we can apply Proposition 13.11 to obtain a sequence, {p an }n C {p a }a>o, with a„ — > 0, and a Vishik- 
Fursikov measure p over / such that p an p in Vtyi). First, let us check that p is an invariant Vishik-Fursikov 
measure, which by ([5]) is equivalent to show that 

/ ^(cr T u)dp(u) = / 93(u)dp(u), for all (p G C b {yi). 
Jyi Jyi 

We already have that for all <p <E Cb(yi), 

/ tp(a T u)dp an (u) = <^(u)dp Q „(u). 

Jyi Jyi 

Note also that <p o a T belongs to Cbiyi) for all <p & Cb(yi), so that, since p an — p in V(yi), we obtain that 

/ <^(cr T u)dp(u) = lim / <^(cr T u)dp Q „(u) = lim / ip(u)dp an (u) = / ^(u)dp(u). 

Jy, ^ooJ yi n^e C Jy i J yi 

Define p := TltP, for an arbitrary time t £ I. Then, p is a stationary Vishik-Fursikov statistical solution by 
definition. Moreover, p an — ^ \i in V(H W ). Indeed, let € Cb(H w ), then y> := o Tl t belongs to Cbtyi) and 

lim / cj)(u)dp an (u) = lim / (p(u)dp an (u) = / <^(u)dp(u) = / 0(u)d/x(u). 

n^ooJ H n-yooj yi Jy { J H 

This concludes the proof. □ 
4 Conclusions 

We have proved that, under natural conditions, families of statistical solutions of the 3D Navier-Stokes-a model, 
depending on the parameter a, possess, as a goes to zero, subsequences that converge to a statistical solution of 
the 3D Navier-Stokes equations. The main condition for the convergence is that the mean kinetic energy of the 
family of statistical solutions be uniformly bounded. This yields suitable tightness and compactness properties 
needed for the existence of a convergent subsequence. 

The statistical solutions contain the statistical information of a given flow, including turbulent flows, which 
are the case of most interest. It is therefore natural to study conditions that guarantee that the statistical 
information obtained from an approximate problem converges to the original problem in a suitable sense. 
We succeeded in proving this convergence under simple and natural conditions. This result implies that the 
statistical information obtained from the 3D Navier-Stokes-a model are good approximations of the statistical 
information of flows modelled by the 3D Navier-Stokes equations. 

Moreover, the techniques developed in this paper seem to allow for an extension of the result to a wide range 
of models and approximations. This is currently a work in progress and will be presented elsewhere. 
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